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(57) Abstract: A pre-distorter pre-distorts an input signal prior to being applied to an amplifier in order to reduce spurious emissions 
in the resulting amplified signal. The pre-distorter implements an inverted version of a model of the amplifier that models both the 

^ frequency -independent (Fl) characteristics of the amplifier as well as the frequency-dependent (FD) characteristics of the amplifier. 

^ Techniques and architectures are presented for (1) generating and updating the model, (2) inverting the model, and (3) updating the 

^* inverted model. 
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ADAPTIVE DIGITAL PRE-DISTORTION USING AMPLIFIER MODEL 
THAT INCORPORATES FREQUENCY-DEPENDENT NON-LINEARITIES 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to signal processing, and, in particular, to the pre-distortion of 
signals for transmission, for example, in a wireless communication network, to reduce spurious 
emissions. 

Cross-Reference to Related Applications 

This application claims the benefit of the filing date of provisional application no. 60/369,488, 
filed on 04/03/02 as attorney docket no. C0013PROV. 

Description of the Related Art 

Modem wireless communication networks employ complex modulation schemes that necessitate 
tight control of spurious emissions (sometimes called "out-of-band emissions") in order to avoid 
interfering with neighboring carriers and to comply with the requirements of regulatory bodies (e.g., 
FCC) and standards bodies (e.g. ITU). One source of spurious emissions is the base station transmitter 
amplifier that is used to amplify signals prior to transmission as wireless (e.g., RF) signals to wireless 
(e.g., mobile) units in a wireless communication network, such as a cellular voice and/or data network. 
Prior art techniques for reducing such spurious emissions were able to satisfy previous requirements. 
However, recent developments in wireless communication networks (e.g., Universal Mobile 
Telecommunication Service (UMTS)) place additional burden on the base station transmitter amplifier 
and make it advantageous to reduce the spurious emissions even further. 

In the prior art, quasi-linear amplifiers operating in Class A, A/B, or B modes have been modeled 
as a memory-less non-linearity. (As is known in the art, different classes of operation relate to different 
quiescent operating points of the amplifiers. In Class A, the transistor is "on" for the whole cycle of a 
sine wave input. In Class B, the transistor is "on" for half the cycle. And, in Class A/B, the transistor is 
one for more than 50%, but less than 100% of the cycle.) Under the assumption of being modeled as a 
memory-less non-linearity, the input-output relationship is given by Equation (1) as follows: 

y = G{a x )x 

= (G^H JG g (a x ))x (i) 

= G m {a x )e jG ' M x 
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where: 

X = d x €^ x is the complex baseband representation of the input signal, where a x is the 

amplitude and <j> x is the phase of the input signal x , 

y is the complex baseband representation of the output signal, and 

G/(0» G q (*) are arbitrary (i.e., unspecified but suitable) functions of the input signal envelope. 

In the above model, it is assumed that the instantaneous gain of the amplifier is solely a function of the 
instantaneous input envelope. This model can be used to describe both AM/AM and AM/PM distortions 
that are observed in most amplifiers operating in Class A, A/B, and B modes, where AM means 
amplitude modulation and PM means phase modulation. A digital pre-distorter can be constructed to 
linearize an amplifier that is described by the above model. 

Fig. 1 shows a block diagram of a conventional amplifier pre-distortion architecture 100 
comprising digital pre-distorter 102, which pre-distorts the input signal prior to amplification by 
memory-less amplifier 104. Pre-distorter 102 implements: 

1 . A method for computing the instantaneous input envelope a x , and 

2. A method of pre-distorting the input signal by multiplying the input signal by a complex gain that 
is solely a function of the instantaneous input envelope. 

If the pre-distorted signal is given by X = G pd (a x )x where G pd (•) is the pre-distorter gain, then the 

input-output relationship for the cascaded pre-distorter and amplifier system can be written according to 
Equation (2) as follows: 

y = G{a x )x 

= Gi\G pi {a x )\a x )G p M* ® 

The pre-distorter gain is computed so that the gain of the pre-distorter cascaded with the amplifier is 
constant. Therefore, for the ideal pre-distorter G^G pd (a x )\a x ^G pd {a x ) = G tgt where G tgt is the 

target gain for the amplifier. 

Fig. 2 is a block diagram of amplifier pre-distortion architecture 100 of Fig. 1, showing further 
detail on a conventional implementation of pre-distorter 102 in which pre-distorter 102 pre-distorts the 
input signal prior to amplification by memory-less amplifier 104. In particular, pre-distorter 102 
comprises a model 202 of amplifier 104. The input signal x is applied to the amplifier model 202 to 



WO 03/085822 PCT/US03/07046 

-3- 

generate a model G(a x )x of the distorted output signal. Difference node 204 generates an estimate of 

the input-output error based on the input signal X and the modeled distorted output signal G(a x )x 
according to Equation (3) as follows: 

e^x) = G(a x )x-x, (3) 

where G(-) is an estimate of the amplifier gain. The pre-distorted signal X is generated at difference 
node 206 by subtracting the estimated error e~ x from the input signal X , where the pre-distorted signal 

X is then applied to amplifier 104. 

For the sake of clarity, assume that the input and output signals are normalized such that the 
target gain is unity. Further, assume that the amplifier behaves as a true memory-less linearity. It can be 
shown that, if amplifier model 202 is constructed accurately, the output y of amplifier 104 with the 

pre-distorted signal input X is given by Equation (4) as follows: 

y = G(a x )x 
* G(a x )x-e 9x 

(4) 



= x+x{G(a x )-G{a x j) 



« x 



Memory-less amplifier model 202 is only an approximate model for most amplifier systems. 
Consequently, the optimal pre-distorter constructed from such a model cannot completely linearize most 
amplifiers. 

SUMMARY OF THE INVENTION 
The problems in the prior art are addressed in accordance with the principles ofthe present 
invention by an amplifier model that is able to model the performance of most Class A, A/B, and B 
amplifiers more accurately than the simple memory-less models ofthe prior art, such as model 202 of 
Fig. 2. In certain embodiments of the present invention, a pre-distorter pre-distorts an input signal prior 
to being applied to an amplifier in order to reduce spurious emissions in the resulting amplified signal. 
The pre-distorter implements an inverted version of a model of the amplifier that models both the 
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frequency-independent (FI) characteristics of the amplifier as well as the frequency-dependent (FD) 
characteristics of the amplifier. 

In one embodiment, the present invention is a method and apparatus for pre-distorting a signal 
for amplification. According to the embodiment, an input signal is received, and pre-distortion is applied 
to the input signal to generate a pre-distorted signal, such that, when the pre-distorted signal is applied to 
an amplifier to generate an amplified signal, the pre-distortion reduces spurious emissions in the 
amplified signal. The pre-distortion is generated using an inverse of a model of the amplifier, where the 
model comprises a model of frequency-independent characteristics of the amplifier in combination with a 
model of frequency-dependent characteristics of the amplifier. 

In another embodiment, the present invention is an apparatus for pre-distorting a signal for 
amplification. Each of one or more high-order transfer function elements is configured to apply a 
transfer function of a different order greater than one to a corresponding distortion product for an input 
signal. A summation node is configured to sum the input signal and the output of each high-order 
transfer function element. An inverted transfer function element is configured to apply an inverted first- 
order transfer function to the output of the summation node. An inverted FI element is configured to 
invert the frequency-independent gain of the amplifier to generate a pre-distorted signal, such that, when 
the pre-distorted signal is applied to an amplifier to generate an amplified signal, spurious emissions in 
the amplified signal are reduced. 

In another embodiment, the present invention is a method for generating a model of an amplifier, 
wherein the model comprises an FI model of frequency-independent characteristics of the amplifier in 
combination with an FD model of frequency-dependent characteristics of the amplifier, wherein the 
model is generated by estimating a plurality of transfer functions for the FD model in a sequence 
corresponding to increasing order of the transfer function. 

In another embodiment, the present invention is a method for adaptively updating a pre-distorter 
configured to pre-distort an input signal for application to an amplifier configured to generate an 
amplified signal, the pre-distorter comprising at least one look-up table (LUT), wherein the at least one 
LUT is adaptively updated by minimizing an error between the input signal and a feedback signal based 
on the amplified signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other aspects, features, and advantages of the present invention will become more fully apparent 
from the following detailed description, the appended claims, and the accompanying drawings in which 
like reference numerals identify similar or identical elements. 

Fig. 1 shows a block diagram of a conventional amplifier pre-distortion architecture; 

Fig. 2 is a block diagram of the amplifier pre-distortion architecture of Fig. 1, showing further 
detail on a conventional implementation of the pre-distorter; 



WO 03/085822 PCT/US03/07046 

-5- 

Fig. 3 shows a block diagram of an amplifier pre-distortion architecture, according to one 
embodiment of the present invention; 

Fig. 4 shows a discrete-time representation of the model of a Class A, A/B, or B amplifier, 
according to one embodiment of the present invention; 

Fig. 5 shows a discrete-time representation of an amplifier model, according to an alternative 
implementation of the present invention; 

Fig. 6 shows a flow diagram of the processing implemented by the pre-distorter of Fig. 3, 
according to one embodiment of the present invention; 

Fig. 7 shows an architecture that can be used to invert the amplifier model of Fig. 4, according to 
one embodiment of the present invention; 

Fig. 8 shows a high-level block diagram illustrating the operation of the pre-distorter of Fig. 3; 

and 

Fig. 9 shows a block diagram of an architecture for the pre-distorter of Fig. 3, according to one 
embodiment of the present invention. 

DETAILED DESCRIPTION 
Fig. 3 shows a block diagram of an amplifier pre-distortion architecture 300, according to one 
embodiment of the present invention. As shown in Fig. 3, the discrete-time, input signal x[n] is applied 

to digital pre-distorter 302, which generates a digital, pre-distorted signal p(x[n]) based on a model of 
high-power amplifier (HPA) 304. The digital, pre-distorted signal then passes through channel 303 and 
arrives at amplifier 304 as an analog, pre-distorted signal p(x(t)) . Channel 303 includes appropriate 
digital-to-analog converters, IF to RF upconverters, mixers, and filters that convert the digital, pre- 
distorted signal p(x[n]) into the analog, pre-distorted signal p(x(t)) . Amplifier 304 generates 
analog, output signal h(p(t)) , where A() is the transfer function of amplifier 304 whose gain is 
K G . Architecture 300 also includes a feedback (FB) channel 305, which samples the output signal 

h(p(t)) and provides a digital FB signal h(p[n]) to pre-distorter 302, which periodically updates its 
processing according to that feedback signal. 

Fig. 4 shows a discrete-time representation of model 400 of a Class A, A/B, or B amplifier, such 
as amplifier 304 of Fig. 3, according to one implementation of the present invention. As indicated in Fig. 
4, the output y[n] of the amplifier is modeled as a sum of several components y[n] , where: 



WO 03/085822 PCT/US03/07046 

-6- 

o y x [n] is the result of applying the input signal x[n] to model 402 and then to first-order 
transfer function H x (z) module 404-1, which models the output matching networks. Model 
402, which models AM/AM and AM/PM distortion to the first order, applies a gain G(a x ) , 
which is a memory-less nonlinear function of the input signal envelope a x ; 
5 o j> 2 [n ] is the output of an arbitrary transfer function acting on a second-order distortion product 

X 2 [n] of the input according to Equation (5) as follows: 

y 2 [n]^al[n]e JMtt] *h 2 [n] (5) 

where " * " represents convolution and h 2 [n] is the impulse response corresponding to second- 
order transform function H 2 (z) module 404-2, which applies the z-transform of the transfer 

10 function; 

o y 3 [n] is the output of an arbitrary transfer function acting on a third-order distortion product 

JC 3 [n] of the input according to Equation (6) as follows: 

y,[n] = al[n]e j * An] *h,[n] (6) 

where [n] is the impulse response corresponding to third-order transform function H 3 (z) 

1 5 module 404-3, which applies the z-transform of the transfer function; and so on until 

o y N [n] is the output of an arbitrary transfer function acting on a JV^-order distortion product 

X N [ri] of the input according to Equation (7) as follows: 

yM = a?[n)e»< [ " ] *h N [n] (7) 

where h N [n] is the impulse response corresponding to A^-order transform function H N (z) 

20 module 404-N, which applies the z-transform of the transfer function. 

These W components are added together at summation node 406 to form the modeled amplifier 

output signal y[n] , where, in general, N may be any positive integer. Model 400 becomes more 
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accurate as more higher-order products are added to the output. Note that the transfer function H } (z) , 
whose impulse response is denoted by > models the frequency dependence resulting from the 

output matching networks and combiners. For most, practical applications, the transfer functions 
H i (z) can be modeled as causal finite impulse response (FIR) filters. 

5 Fig. 5 shows a discrete-time representation of an amplifier model 500, according to an alternative 

implementation of the present invention. In model 500, AM/AM-AM/PM model 502 and first-order 
transfer function module 504-1 are analogous to AM/AM-AM/PM model 402 and first-order transfer 
function module 404-1 of Fig. 4. In model 500, however, all the higher-order distortion products of the 
input signal are processed by a single, second-order transfer function H 2 (z) module 504-2. 

10 Furthermore, as shown in Fig. 5, the input signal x[n] is subtracted from the pre-distorted 

output of AM/AM-AM/PM model 502 at difference node 503. As such, the output y 2 [n] from second- 
order transfer function module 504-2 may be given according to Equation (8) as follows: 

y 2 [n] = (G (\x[n]\ - x[n]) x[n) * \ [«]) (8) 

Amplifier model 500 should be adequate for most pre-distortion applications. If better performance is 
1 5 required from the pre-distorter, the more general amplifier model 400 of Fig. 4 can always be used. 

Independent of which amplifier model of the present invention is to be used, the parameters for 
the selected model need to be computed. This includes computing estimates of the memory-less 
amplifier gain GQ and all of the transfer functions H ( (z) . Once computed, an amplifier model in 

accordance with the present invention can be substituted for prior-art amplifier model 202 in the pre- 
20 distortion architecture of Fig. 2, and the resulting pre-distorter architecture (or any other suitable 
architecture) can be used to invert the model and linearize the amplifier. 

Fig. 6 shows a flow diagram of the processing involved in generating and updating an amplifier 
model, such as model 400 of Fig. 4 or model 500 of Fig. 5, according to one embodiment of the present 
invention. At step 602 of Fig. 6, the gain of the amplifier is estimated. If actual time-synchronized 

25 samples of the input ( X = [ JC, JC 2 • • • X N ]' ) and output ( y = [y x y 2 • • • y N ]' ) signals at baseband are 



available, then the measured gain G red can be calculated according to Equation (9) as follows: 
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'measured 



x x 2 



(9) 



where [•]' represents the transpose of a column vector. An arbitrary linear function 0(1*1) can then 

be fit to the measured data to model the amplifier gain. The parameters for the arbitrary function can be 
estimated by the method of least squares. That is, the parameters of the function that models the gain are 
computed by minimizing the cost function % Si ven Equation (10) as follows: 

l|2 



r measured 



-G 



(10) 



IIW 



where G = [g(|xJ) G(|a: 2 |) G(|jc^|)]' and W is a weighting matrix. The choice of the 

weighting matrix may depend on the specific application. For example, the weighting matrix can be 
chosen to weight the error differently at the higher power levels than at the lower power levels. If all 
10 power levels are to be weighted equally, then the weighting matrix can be chosen as the identity matrix. 
Note that any arbitrary nonlinear function can be fit to the measured data. Polynomials and splines are 
special cases of arbitrary nonlinear functions that can be used. The choice of the appropriate nonlinear 
function should be apparent to those skilled in the art. 

After the amplifier gain has been estimated in step 602, the transfer functions are estimated in 

1 5 step 604. In particular, the transfer function H x (z) may be estimated by minimizing the cost function 
% 0 given by Equation (1 1) as follows: 

% 0 =E{y[n]- /*,[«]* x,[n]) 2 (11) 

where x x [n] = G(a x [n])x[n] is the output of the memory-less AM/AM-AM/PM model (e.g., model 

402 of Fig. 4 or model 502 of Fig. 5). The optimum solution can be obtained by employing any of the 
20 known adaptive-estimation techniques, such as recursive least-squares (RLS), least mean square (LMS), 
or least-squares solutions. 

After estimating H x (z) , the transfer function H 2 (z) may be estimated by minimizing the cost 

function^ given by Equation (12) as follows: 
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Xx =^W»]-AW-*iW*^W) a - (12) 

Once again, the optimum solution can be obtained by employing any of the known adaptive-estimation 
techniques. Alternatively, the solution can be obtained in the frequency domain as follows. The 
optimum solution for H 2 (z) at a given frequency can be obtained by minimizing the cost function in 
5 the frequency domain and weighting the frequency points of interest. This can then be repeated for 
several frequency points. Once the frequency response of H 2 (z) has been estimated at N distinct 

frequency points, the time-domain impulse response { h 2 [n]} can be obtained by a least-squares fit to 

the desired frequency response. 

Next, the transfer function H 3 (z) may be estimated by minimizing the cost function Xi & yeYi 

10 by Equation (13) as follows: 

X 2 ^E(y[n]-Un]-^[n]*x 3 [n]f 03) 

where j> 2 [rt] = j>, [n] + h 2 [n] * X 2 [n] . Any higher-order transfer functions H i (z) , for i > 3 , can 

be computed in analogous fashion. 

Once the transfer functions have been computed, an inverted version of the amplifier model can 
15 be constructed and used to generate the pre-distorted signal (step 606) that is then applied to the 

amplifier (step 608). The pre-distorter architecture shown in Fig. 2 (or any other suitable architecture 
that inverts the amplifier model of Fig. 4 or Fig. 5) can be used to pre-distort the signal and linearize the 
amplifier. 

Fig. 7 shows an architecture 700 that can be used to invert amplifier model 400 of Fig. 4, 
20 according to one embodiment of the present invention. In practice, FIR filters can be used to 

approximate the transfer functions in architecture 700. As suggested in Fig. 7, architecture 700 may be 
implemented with any one or more paths of negated transfer functions of order two or greater. It is also 
possible to implement architecture 700 without any negated transfer function paths. In that case, the 
summation node may be omitted, where the input signal alone is input to the inverted first-order transfer 
25 function, followed by the inverse of the AM/AM-AM/PM model. 

The estimates of the amplifier gain and transfer functions generated in steps 602 and 604 can be 
adaptively updated based on the amplifier output signal using suitable recursive estimation techniques. 
This updating is indicated in Fig. 6 by the processing returning from step 608 back to step 602 for 
another pass. This updating can be implemented intermittently, e.g., at fixed intervals or as necessary. 
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Recursive estimation techniques for adaptive filters are well known. See, e.g., Simon Haykin, 
Adaptive Filter Theory, Third Edition, Prentice Hall, 2001, the teachings of which are incorporated 
herein by reference. Many of these techniques can be applied to recursively estimate the transfer 
functions H i (z) . However, recursive estimation techniques for computing the AM/AM-AM/PM pre- 

5 distorter (e.g., corresponding to model 402 of Fig. 4 or model 502 of Fig B) are not that well known. The 
following discussion describes a novel method for recursively computing the memory-less pre-distorter. 

Fig. 8 shows a high-level block diagram illustrating the operation of pre-distorter 302 of Fig. 3, 
according to one implementation of the present invention. In this implementation, pre-distorter 302 
relies on look-up table (LUT) operation presuming the digital input signal x[n] has been converted to 

10 magnitude-phase format. The input magnitude signal is applied in parallel to both AM/AM LUT control 
802 and AM/PM LUT control 804. The output from AM/AM LUT control 802 is the pre-distorted 
magnitude signal, while the pre-distorted phase signal is generated at summation node 806 by summing 
the output from AM/PM LUT control 804 with the input phase signal. 

In preferred implementations, the LUTs for the AM/AM and AM/PM portions of the pre- 

15 distorter operation are adaptively revised on a periodic basis. This revision may be based on minimizing 
the mean-square error between the feedback signal h(p[n]) in Fig. 3 and the input signal x[n] . The 

cost function Xk used in * e feedback control is based on a current observation data set (where the 

input and feedback signals are assumed to have been time-aligned by performing time-delay 
equalization) and is defined according to Equation (14) as follows: 

where S K = (^(1) e(2) • • • €(K)^j , and s(n) is the normalized error between the input sample 
X n and the corresponding feedback amplifier sample (h o p)(x n ) , where (h o p)(x n ) represents 
h [ p[x n ]] . The normalized error is given by Equation (15) as follows: 

£W >„-(*°p)(*j) (15) 
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where h(x) is the amplifier transfer function, p(x) is the pre-distortion polynomial used to generate 

the AM/AM and AM/PM LUTs, and o\ is the sample variance over the observation set 

\x k : k = \,...,K}} . The pre-distortion polynomial p(x) is given by Equation (16) as follows: 

p(x) = f d a m x m =X l * (16) 

where a = (fl 0 , a x 9 ...a N )' is the vector of coefficients for characterizing the LUT generator 
polynomial and X* = (l, X ] , . . . , X N ) is the basis vector for linear space represented by the input data. 

2 

The sample variance G K is given by Equation (17) as follows: 

1 K 2 

^=^—^Z(^-^) O 7 ) 

where m K is the sample mean over the observation set and is given by Equation (18) as follows: 

1 K 

& *=i 

The optimization criteria 2i K for the adaptation are given by Equation (19) as follows: 

a* = argmin{;^} (19) 

a 

where 3Xg 1H1I1 { /(*)} returns the value of x that minimizes the function^*,). The optimization 



2 

criteria 2i K are unique and are obtained by equating the gradient of the cost function Xk res P ect 
to the coefficient vector a to zero as in Equation (20) as follows: 
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where the matrix H is defined by Equation (21) as follows: 

h l 0 - 

1 



H'= — 

0> 



x x n N ■■■ x K n N 



(21) 



where h k m =d m (h° p){x k ) and is defined by Equation (22) as follows: 



d m (h°p)(x) = 



x m 



P(x) 



h • P(x) + r p2(3C) l2 (0? - - p(x)(p 'o />(*)) 

[(/><>/>)(*)] 



(22) 



5 

5 where 3 = — is short-hand notation for a derivative or gradient, depending on the context, and 
2 dz 

p\x) is the first derivative of the polynomial p(x) with respect to X. 

Stochastic approximation can be applied to this problem to successively minimize the gradient of 
Equation (20). This results in the following update procedure for the coefficient set given by Equation 
(23) as follows: 

10 *w=*oU+M*K-* M ) (23) 

where A is a small scalar adaptation gain (e.g., typically about 0.0005) and 

a K = (h'H) 1 H % £ K or * K =T) x s K 

where D is the pseudo-inverse of the matrix H . 

The expression in Equation (22) depends only on the following three things: 
15 (1) The amplifier output feedback observation (h o p)(x K ) = h(p(x K )) corresponding to the 

input signal X K = {x k : k = \ y „. y K} ; 
(2) p(x K ) and p o p(x K ) are obtained from a look-up table and a recursive application of that 
table to itself, respectively; and 
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(3) p 'o p(x K ) which is (or could be) one or more separate LUTs based on the first derivative of 

the polynomial p(x) evaluated on the pre-distorter LUT output 

If the polynomial p(x) is constrained to pass through the origin a Q = 0 , then the top row in the 

matrix H can be deleted. The rest of this derivation relies on that assumption. 
5 Although the SVD is the preferred method of solution from the standpoint of numerical stability, 

it has the following three main drawbacks: 

(1) The computation cost is significant, being proportional to K 2 ; 

(2) Storage requirements are larger than other methods; and 

(3) In general, double-precision computation should be used for many of the intermediate steps. 

10 As an alternative, the pseudo-inverse of H may be calculated directly using 0(N 2 + KN) FLOPS, 

provided the matrix is not allowed to become ill-conditioned. One approach, which is effective against 
this condition, is based on the following processing: 
Define a scalar y by the algorithm: 

if{z 2 K>zL) 

15 theny = 5y 

else y = 0.8^ 
endif 

2 2 

The value for % last is the previous iterations mean square error % K . 
The product matrix H*H is then (re-)defined as 

20 H'H <r- H'H + X <flag(H'H) 

Using this revised definition, the pseudo-inverse is calculated by inverting the NxN matrix H*H . 
By applying this procedure, a stable and more-efficient method may be implemented that should provide 
a 'faster* stable solution to the update procedure. 

The final step in the update procedure is to recalculate and replace the AM/AM and AM/PM 

25 LUTs (and the "derivative" table) using the revised set of polynomial coefficients a new , and to replace 

the previous coefficient set with this newer estimate, 2i old = a ncw . The derivative LUT generator 
polynomial is given by Equation (24) as follows: 
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P'(x)=S^ M - , w 

where the coefficients are the same as those for the pre-distortion LUTs. 

In other implementations, it is possible that the pre-distorter is implemented with only one LUT 
(e.g., either only the AM/AM LUT or only the AM/PM LUT). In such implementations, the above- 
described update method would be applied to only the one LUT. 

Fig. 9 shows a block diagram of architecture 900 for pre-distorter 302 of Fig. 3, according to one 
embodiment of the present invention. In particular, architecture 900 corresponds to a second-order 
implementation of the pre-distorter architecture shown in Fig. 7, where FIR 908 implements the transfer 

function -H? (z) , FIR 912 implements the transfer function , and LUT 918 and multiplier 

920 implement the inverse of the AM/AM-AM/PM model. Delays 902 and 914 are provided to maintain 
synchronization within architecture 900 by accounting for the processing times of their associated 
parallel processing paths. In general, higher-order terms (3 rd and higher) can also be included, if desired, 
using the structure shown in Fig. 7. 

In general, all signals in Fig. 9 are complex (shown as dark lines) except the outputs of 
rectangular-to-polar (R2P) converters 904 and 916. Each R2P converter computes the amplitude of its 
input signal using the well-known CORDIC algorithm. See, e.g., Andraka, R., "A survey of CORDIC 
algorithms for FPGA based computers", Andraka Consulting Group, Inc., and Voider, J., "The CORDIC 
trigonometric computing technique," IRE Trans. Electronic Computing, Vol. EC-8, pp. 330-334, 
September 1 959, the teachings of both of which are incorporated herein by reference. This 
implementation includes the LUT-based correction for the AM/AM-AM/PM effects, the equalizer for the 
channel distortion, and the correction for the second-order, frequency-dependent non-linearity. In 
alternative embodiments, CORDIC R2P converters 904 and 916 could be replaced by other components 
that generate a suitable measure of the input signal, where, depending on the implementation, that 
measure could be the amplitude, power, or even a higher-order measure of the input signal. 

Fig. 9 corresponds to a second-order, baseband implementation of the pre-distorter architecture 
of Fig. 7. Non-baseband (e.g., RF) implementations are also possible. In one possible second-order, RF 
implementation, FIR 912 is replaced by a baseband-to-RF converter followed by an RF filter, and 
multiplier 920 is replaced by a vector modulator. 

The present invention may be implemented in the context of wireless signals transmitted from a 
base station to one or more mobile units of a wireless communication network. In theory, embodiments 
of the present invention could be implemented for wireless signals transmitted from a mobile unit to one 
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or more base stations. The present invention can also be implemented in the context of other wireless 
and even wired communication networks to reduce spurious emissions. 

Embodiments of the present invention may be implemented as circuit-based processes, including 
possible implementation on a single integrated circuit. As would be apparent to one skilled in the art, 
5 various functions of circuit elements may also be implemented as processing steps in a software program. 
Such software may be employed in, for example, a digital signal processor, micro-controller, or general- 
purpose computer. 

The present invention can be embodied in the form of methods and apparatuses for practicing 
those methods. The present invention can also be embodied in the form of program code embodied in 

10 tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage 
medium, wherein, when the program code is loaded into and executed by a machine, such as a computer, 
the machine becomes an apparatus for practicing the invention. The present invention can also be 
embodied in the form of program code, for example, whether stored in a storage medium, loaded into 
and/or executed by a machine, or transmitted over some transmission medium or carrier, such as over 

15 electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the 
program code is loaded into and executed by a machine, such as a computer, the machine becomes an 
apparatus for practicing the invention. When implemented on a general-purpose processor, the program 
code segments combine with the processor to provide a unique device that operates analogously to 
specific logic circuits. 

20 It will be further understood that various changes in the details, materials, and arrangements of 

the parts which have been described and illustrated in order to explain the nature of this invention may be 
made by those skilled in the art without departing from the scope of the invention as expressed in the 
following claims. 
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CLAIMS 

What is claimed is: 

1 . A method for pre-distorting a signal for amplification, the method comprising the steps of: 
(a) receiving an input signal; and 

5 (b) applying pre-distortion to the input signal to generate a pre-distorted signal, such that, when the 

pre-distorted signal is applied to an amplifier to generate an amplified signal, the pre-distortion reduces 
spurious emissions in the amplified signal, wherein: 

the pre-distortion is generated using an inverse of a model of the amplifier; and 
the model comprises a model of frequency-independent (FI) characteristics of the amplifier (an FI 
10 model) in combination with a model of frequency-dependent (FD) characteristics of the amplifier (an FD 
model). 

2. The invention of claim 1 , further comprising the step of generating at least one of the FI model 
and the FD model. 

15 

3 . The invention of claim 2, wherein the FD model is generated by estimating a plurality of transfer 
functions in a sequence corresponding to increasing order of the transfer function. 

4. The invention of claim 2, wherein the model is generated by: 

20 (1) estimating the FI model based on an estimate of the frequency-independent gain of the amplifier; 

(2) estimating a first-order transfer function for the FD model using the estimated FI model; and 

(3) estimating a second-order transfer function for the FD model using the estimated FI model and 
the estimated first-order transfer function. 

25 5. The invention of claim 4, wherein: 

the FI model is estimated by minimizing a cost function that is dependent on the gain of the 
amplifier; 

the first-order transfer function is estimated by minimizing a cost function that is dependent on the 
estimated FI model; and 

30 the second-order transfer function is estimated by minimizing a cost function that is dependent on the 

estimated FI model and the estimated first-order transfer function. 



35 



6. The invention of claim 4, wherein the model is further generated by estimating one or more 
transfer functions for the FD model of order greater than two, wherein a transfer function of order n is 
estimated using the estimated FI model and each estimated transfer function of order less than n. 
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7. The invention of claim 2, further comprising the step of adaptively updating the model. 

8. The invention of claim 7, wherein the model is adaptively updated by adaptively updating a 
plurality of transfer functions in a sequence corresponding to increasing order of the transfer function. 

9. (e.g., Figs. 4 & 5) The invention of claim 1 , wherein the FD model comprises a plurality of 
different-order transfer functions, each of which is configured to operate on a different-order distortion 
product based on the input signal, wherein: 

the output from the FI model is applied to a first-order transfer function; and 
the outputs from the plurality of transfer functions are summed to generate an estimate of the 
amplified signal. 

10. (e.g., Fig. 5) The invention of claim 9, wherein the FD model comprises a second-order transfer 
function that is configured to operate on a difference between the input signal and the output from the FI 
model. 

1 1 . (e.g., Fig. 4) The invention of claim 9, wherein the FD model comprises one or more transfer 
functions having order two or greater, each of which is configured to operate on a distortion product 
having a corresponding order of two or greater. 

12. (e.g., Fig. 4) The invention of claim 1 1, wherein the FD model comprises two or more transfer 
functions having order two or greater, each of which is configured to operate on a distortion product 
having a corresponding order of two or greater. 

13. The invention of claim 9, wherein, in generating the model, the plurality of transfer functions are 
estimated in a sequence corresponding to increasing order. 

14. (e.g., Fig. 7) The invention of claim 1, wherein the inverse of the model comprises: 

(1) zero, one, or more negated high-order transfer function elements, each corresponding to a 
negation of one of one or more transfer functions of the FD model of order greater than one and 
configured to operate on a corresponding distortion product; 

(2) an inverted first-order transfer function element corresponding to an inverse of a first-order 
transfer function of the FD model and configured to operate on the sum of the input signal and the zero, 
one, or more outputs from the zero, one, or more negated high-order transfer function elements; and 

(3) an inverted FI model element corresponding to an inverse of the FI model and configured to 
operate on the output from the inverted first-order transfer function to generate the pre-distorted signal. 
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15. (e.g., Fig. 9) The invention of claim 14, wherein: 

the zero, one or more negated high-order transfer function elements include a negated second-order 
transfer function element comprising: 

(A) a first component configured to generate a measure of the input signal; 

(B) a first multiplier configured to multiply the input signal by the measure of the input signal to 
generate a second-order distortion product; and 

(C) a first filter configured to apply a negated second-order transfer function to the second-order 
distortion product; 

the inverted first-order transfer function element comprises a second filter configured to apply the 
inverted first-order transfer function to the sum of the input signal and the one or more outputs from the 
one or more negated high-order transfer function elements; and 

the inverted FI model element comprises: 

(A) a second component configured to generate a measure of the output from the second filter; 

(B) a look-up table configured to retrieve one or more predistortion components based on the 
measure of the output from the second filter; and 

(C) a second multiplier configured to multiply the output from the second filter by the one or 
more predistortion components to generate the pre-distorted signal. 

16. The invention of claim 15, wherein: 

the first component is configured to generate a measure of the amplitude or the power of the input 
signal; and 

the first filter is a finite impulse response (FIR) filter. 

17. The invention of claim 16, wherein the second filter is an FIR filter corresponding to the inverse 
of a first-order transfer function. 

18. The invention of claim 16, wherein: 

the inverted first-order transfer function further comprises a baseband-to-RF converter preceding the 
second filter; 

the second filter is an RF filter corresponding to the inverse of a first-order transfer function; and 
the second multiplier is a vector modulator. 

19. (e.g., Fig. 8) The invention of claim 1, wherein: 

the inverse of the model comprises an inverse of the FI model; 

the inverse of the FI model generates an inverted FI model output having an amplitude and a phase; 
the magnitude of the inverted FI model output is a function of the magnitude of the input signal; and 
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the phase of the inverted FI model output is a function of both the magnitude and phase of the input 
signal. 

20. The invention of claim 19, wherein: 
5 the magnitude of the inverted FI model output is derived using the input signal magnitude as an index 

into an AM/AM look-up table; 

the phase of the inverted FI model output is derived by summing the input signal phase and a value 
retrieved using the input signal magnitude as an index into an AM/PM look-up table. 



10 21. The invention of claim 20, wherein at least one of the look-up tables is adaptively updated by 

minimizing an error between the input signal and a feedback signal based on the amplified signal. 

22. The invention of claim 21, wherein the at least one look-up table is updated by applying a 
recursive estimation technique to the look-up table. 

15 

23. The invention of claim 22, wherein the recursive estimation technique is based on non-linear 
least mean squares estimation. 



24. An apparatus for pre-distorting a signal for amplification, wherein the apparatus is configured to: 
20 (a) receive an input signal; and 

(b) apply pre-distortion to the input signal to generate a pre-distorted signal, such that, when the pre- 
distorted signal is applied to an amplifier to generate an amplified signal, the pre-distortion reduces 
spurious emissions in the amplified signal, wherein: 

the pre-distortion is generated using an inverse of a model of the amplifier; and 
25 the model comprises an FI model of frequency-independent characteristics of the amplifier in 

combination with an FD model of frequency-dependent characteristics of the amplifier. 



25. (Fig. 7) An apparatus for pre-distorting a signal for amplification, the apparatus comprising: 

(a) one or more high-order transfer function elements, each configured to apply a transfer function of 
30 a different order greater than one to a corresponding distortion product for an input signal; 

(b) a summation node configured to sum the input signal and the output of each high-order transfer 
function element; 

(c) an inverted transfer function element configured to apply an inverted first-order transfer function 
to the output of the summation node; and 
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(d) an inverted FI element configured to invert the frequency-independent (FI) gain of the amplifier 
to generate a pre-distorted signal, such that, when the pre-distorted signal is applied to an amplifier to 
generate an amplified signal, spurious emissions in the amplified signal are reduced. 

26. (e.g., Fig. 9) The invention of claim 25, wherein: 

the one or more high-order transfer function elements include a second-order transfer function 
element comprising: 

(A) a first component configured to generate a measure of the input signal; 

(B) a first multiplier configured to multiply the input signal by the measure of the input signal to 
generate a second-order distortion product; and 

(C) a first filter configured to apply a second-order transfer function to the second-order 
distortion product; 

: the inverted first-order transfer function element comprises a second filter configured to apply the 
inverted first-order transfer function to the sum of the input signal and the outputs from the one or more 
high-order transfer function elements; and 
the inverted FI element comprises: 

(A) a second component configured to generate a measure of the output from the second filter; 

(B) a look-up table configured to retrieve one or more predistortion components based on the 
measure of the output from the second filter; and 

(C) a second multiplier configured to multiply the output from the second filter by the one or 
more predistortion components to generate the pre-distorted signal. 

27. The invention of claim 26, wherein: 

the first component is configured to generate a measure of the amplitude or the power of the input 
signal; and 

the first filter is a finite impulse response (FIR) filter. 

28. The invention of claim 27, wherein the second filter is an FIR filter corresponding to the inverse 
of a first-order transfer function. 

29. The invention of claim 27, wherein: 

the inverted first-order transfer function further comprises a baseband-to-RF converter preceding the 
second filter; 

the second filter is an RF filter corresponding to the inverse of a first-order transfer function; and , 
the second multiplier is a vector modulator. 
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30. A method for generating a model of an amplifier, wherein the model comprises an FI model of 
frequency-independent characteristics of the amplifier in combination with an FD model of frequency- 
dependent characteristics of the amplifier, wherein the model is generated by estimating a plurality of 
transfer functions for the FD model in a sequence corresponding to increasing order of the transfer 

5 function. 

3 1 . The invention of claim 30, wherein the model is generated by: 

(1) estimating the FI model based on an estimate of the frequency-independent gain of the amplifier; 

(2) estimating a first-order transfer function for the FD model using the estimated FI model; and 
10 (3) estimating a second-order transfer function for the FD model using the estimated FI model and 

the estimated first-order transfer function. 

32. The invention of claim 31, wherein: 

the FI model is estimated by minimizing a cost function that is dependent on the gain of the 
15 amplifier; 

the first-order transfer function is estimated by minimizing a cost function that is dependent on the 
estimated FI model; and 

the second-order transfer function is estimated by minimizing a cost function that is dependent on the 
estimated FI model and the estimated first-order transfer function. 

20 

33. The invention of claim 31, wherein the model is further generated by estimating one or more 
transfer functions for the FD model of order greater than two, wherein a transfer function of order n is 
estimated using the estimated FI model and each estimated transfer function of order less than n. 

25 34. The invention of claim 30, further comprising the step of adaptively updating the model. 

35. The invention of claim 34, wherein the model is adaptively updated by adaptively updating a 
plurality of transfer functions in a sequence corresponding to increasing order of the transfer function. 

30 36. A method for adaptively updating a pre-distorter configured to pre-distort an input signal for 

application to an amplifier configured to generate an amplified signal, the pre-distorter comprising at 
least one look-up table (LUT), wherein the at least one LUT is adaptively updated by minimizing an 
error between the input signal and a feedback signal based on the amplified signal. 



35 37. The invention of claim 36, wherein the at least one look-up table is updated by applying a 

recursive estimation technique to the at least one look-up table. 
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38. The invention of claim 37, wherein the recursive estimation technique is based on non-linear 
least mean squares estimation. 
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